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Abstract 

We derive Euler-Lagrange type equations for fractional action-like inte- 
grals of the calculus of variations which depend on the Riemann-Liouville 
derivatives of order (a, j3), a > 0, (3 > 0, recently introduced by J. Cresson 
and S. Darses. Some interesting consequences are obtained and discussed. 
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1 Introduction 

Fractional derivatives and integrals play a leading role in the understanding 
of complex, classical or quantum, conservative or nonconservative, dynamical 
systems with holonomic as well as with nonholonomic constraints. Many appli- 
cations of fractional differentiation and integration can be found in turbulence 
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and fluid dynamics, chaotic dynamics, solid state physics, chemistry, stochastic 
dynamical systems, plasma physics and controlled thermonuclear fusion, kinetic 
theory, quantization, field theory, nonlinear control theory, image processing, 
nonlinear biological systems, material sciences, rheological properties of rocks, 
scaling phenomena, astrophysics, etc. (see e.g. [3J HHJ HI]). Its origin goes 
back more than three centuries, when in 1695 L'Hopital made some remarks 
to Leibniz about the mathematical meaning of a fractional derivative of order 
1/2. Leibniz's response was "an apparent paradox, from which one day use- 
ful consequences will be drawn" . In these words fractional calculus was born. 
It was primarily a study reserved to the best minds in mathematics, includ- 
ing J. Fourier, N. H. Abel, J. Liouville, B. Riemann, H. Holmgren, etc., who 
contributed strongly to the fractional analysis program. A deep research was 
carried out by Liouville from 1832 to 1837, where he succeeded in defining the 
first fractional integration operator. Later on, further developments lead to the 
construction of the well-known Riemann-Liouville fractional integral operator, 
which plays nowadays a leading and important role in the analysis of complex 
dynamical systems. There exist today many different forms of fractional integral 
operators, ranging from divided-difference types to infinite-sum types, including 
the Grunwald-Letnikov fractional derivative, the Caputo fractional derivative, 
etc. One can say, however, that the Riemann-Liouville operator is still one of 
the most frequently used for fractional integration. 

The fractional calculus was considered a theoretical mathematical field with 
no physical applications for more than three centuries. Last decades have shown 
the contrary. Several fields of application of fractional differentiation and frac- 
tional integration are already well established, some others have just started. 
One of the areas that recently emerged within the fractional framework and 
which is being subject to strong research is the Calculus of Variations (CoV) 
and respective Euler-Lagrange type equations. In 1996 F. Riewe formulated 
the CoV problem with fractional derivatives and obtained the respective Euler- 
Lagrange equations (ELe's) , combining both conservative and non-conservative 
cases [T7]. In 2001 another approach was developed by M. Klimek by consider- 
ing fractional problems of the CoV but with symmetric fractional derivatives. 
Correspondent ELe's were obtained, using both Lagrangian and Hamiltonian 
formalisms [12]. In 2002 O. Agrawal extended Klimek variational problems by 
considering right and left fractional derivatives in the Riemann-Liouville sense 
[1]. In 2004 Agrawal ELe's were used by D. Baleanu and T. Avkar to investigate 
problems with Lagrangians which are linear in the velocities [2] . In all the above 
mentioned studies, ELe's depend on left and right fractional derivatives, even 
when the problem depends only on one type of them. To avoid this situation, in 
2005 M. Klimek studied problems depending on symmetric derivatives, proving 
that in this case ELe's include only the derivatives that appear in the formu- 
lation of the problem [T3]. One major problem with all the above mentioned 
approaches is the presence of a non-local fractional differential operator with an 
adjoint which is not its symmetric. Other difficulties arise due to a complicated 
Leibniz rule and the non-existence of a fractional analogue to the chain rule. 

The physical reasons for the appearance of fractional equations are, in gen- 
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eral, long-range dissipation and non-conservatism. For this reason, it seems 
of interest to study the fractional Hamiltonian of nonconservative dynamical 
systems. With this in mind, the first author has proposed in 2005 a novel 
approach, entitled Fractional Action-Like Variational Approach (FALVA), to 
model nonconservative dynamical systems [BJ 13 • The Euler-Lagrange equations 
proved to be similar to the classical ones, with no fractional derivatives appear- 
ing, but with the presence of a fractional generalized external force acting on 
the system. The momentum, the Hamiltonian and Hamilton's equations were 
proved to depend on the fractional order of integration and to vary inversely to 
time. In the present work we extend FALVA formalism and previous results by 
considering Riemann-Liouville fractional derivatives of order (a,/3). Respective 
Euler-Lagrange equations are obtained using a Taylor-Riemann series expan- 
sion; the concept of fractional constant of motion is introduced and illustrated. 

2 Fractional Euler-Lagrange Equations 

We begin by recalling the standard definitions of Riemann-Liouville fractional 
integrals and fractional derivatives. 

Definition 2.1 (Left Riemann-Liouville fractional integral). If f {t) & C (a, 6), 
and a > 0, then the left Riemann-Liouville fractional integral of order a is 
defined and denoted by 



Definition 2.2 (Right Riemann-Liouville fractional integral). If f (t) ^ C (a, b), 
and a > 0, then the right Riemann-Liouville fractional integral of order a is 
defined and denoted by 



Definition 2.3 (Left and Right Riemann-Liouville fractional derivatives). If 
/ {t) e (a, b), and a > 0, then the left and right Riemann-Liouville fractional 
derivatives of order a, denoted respectively by D"^ and , are defined by 




a 



b 




D^Jit) 



T (n — a) 



1 



n — a—1 



dr 



a 



T {n — a) 



1 



n—a — 1 



dr , 



t 



where n is such that n — 1 < a < n. 
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The classical integrals, and derivatives of order one, are obtained setting 
a = 1. For an introduction to the fractional calculus we refer the reader to 

In this work we consider the fractional derivative operator of order (a, P) 
recently introduced by J. Cresson and S. Darses [S] (see also 

Definition 2.4. Given a, & G M, a < 6, and 7 € C, the fractional derivative 
operator of order (a, /?), a, /3 > 0, is defined by 
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where i = 



Remark 2.1. The operator D'^'^ extends the classical Riemann-Liouville frac- 
tional derivatives: for 7 = — i we have D'^'^ = I?"^; for 7 — i we obtain 

The following Lemma is useful: we make use of ([T]) to prove the fractional 
Euler-Lagrange equations associated to our problem (cf. proof of Theorem 12. 2p . 

Lemma 2.1 ([4J). // /, g e with f{a) ^ f{b) = or g{a) = g{b) = 0, then 

b 

D^'f^f (t) g{t)dt^ - [ f it) Df'j^g (t) dt . (1) 



We are now in conditions to formulate the (a, [3) fractional action-like vari- 
ational problem. 

Definition 2.5. Consider a smooth manifold M and let L be a smooth La- 
grangian function L : x R'' x K — > M, d > 1. For any piecewise smooth path 
q : [a,b] ^ M satisfying fixed boundary conditions q{a) = qa and q(b) = qb we 
define the following fractional action integral: 

b 

S;fa,b) M = ^ / i (^V^ (r) , g (r) , r) [t - rf-' dr . (2) 

a 

The (a, f3) fractional action-like variational problem consists in finding an ad- 
missible q{-) which minimizes ([2]). 

Remark 2.2. We consider two time variables: the intrinsic time r and the ob- 
server time t. This multi-time characteristic is important in applications and is 
the main ingredient of the theory being developed by C. Udriste [21] . 

Remark 2.3. The fractional action integral ([2]) is a generalization of the FALVA 
action integral of O [7l [8] . 

Theorem 12.21 gives a necessary optimality condition for g(-) to be a solution 
of the (a, P) fractional action-like variational problem. 
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Theorem 2.2 ((a,/3) fractional Euler-Lagrange equations). If q : [a,b] M 
is a minimizer of the (a,/3) fractional action-like variational problem (cf. Defi- 
nitionWT^). then 



dL 

\ — a dL 



,D^-^ (r) , g (r) , r) - ^^".^ {D^ 'q (r) , g (r) , r) 



t-T 



(Z?^'^<;(r),g(r),T) (3) 



where D^'".^ represents the fractional derivative with respect to time t. 



Remark 2.4. Theorem 12.21 is an extension of the Euler-Lagrange equations de- 
rived in FALVA: one just needs to choose /3 = 1 in ([3]) to obtain the Euler- 
Lagrange equations of O [7] . 

Proof. We perform a small perturbation of the generalized coordinates as q 
q + eh, £<t:l. As a resuh, Z?^-^ (q + eh) = D^^^q + sD^^'^h and 

b 

a 

which, doing a Taylor expansion of L [D"'^q + eD"'^h, q + eh, r) in e around 
zero, and integrating by parts, imply that 



b 

e 



r(a) 



— {D^-Pq{r),q{r),T) [t-rT-'hir) 



+ it- tT-' -qT (D^^^q (r) D^^l^h (r) , g (r) , r) D'^:^h (r) 



+ {t- rf-' ^ {D^'^q (r) D^^^h (r) , g (r) , r) D'^^^^h (r) 



dr + O (e) 



Making use of Lemma 12.11 and the least action principle we arrive to . □ 

Definition 2.6. A path q : [a, b] AI satisfying equation ([3]) is said to be a 
fractional extremal associated to the Lagrangian L. 

Definition 2.7. The right hand side of 

F^-f^^^{D-^^Mr),q{r),r), (4) 
T — T — t, defines the fractional decaying friction force. 

Remark 2.5. The fractional decaying friction force satisfy the asymptotic prop- 
erty lim = 0. 
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3 The Fractional Hamiltonian Formalism 



We now consider a more general class of fractional optimal control problems: 

b 



'^'"'(Ifc) = j L{u{T),q (r) , r) (t - r) dr > min , 

a 

Z?^'^9(r) = ^(«(r),g(T),r) , 

a, 6 S M, a < 6. In the particular case where {u, q, r) = u (O reduces to 
To obtain a necessary optimality condition to problem ([5]) we introduce the 
augmented action integral (cf. e.g. [S]): 

b 

a 

-p''-^iT)D^-''q{T)\dT (6) 

where p"-'^ is the fractional Lagrange multiplier and the fractional Hamiltonian 
W-/^ is defined by 

Theorem 3.1. If (q,u) is a minimizer of problem ([5]), then there exists a co- 
vector function p'^'P such that the following conditions hold: 



• the fractional Hamiltonian system 



D^J^.^^p-^P {t) = -^ {u (r) , q (r) (r) , r) 



(7) 



-l-rf \' ! — dq 

• the fractional stationary condition 



^^(«(r),g(T),p"''5(T),r) =0. (8) 

Remark 3.1. For the fractional problem of the calculus of variations ([2]) one has 

H'^^^ {u, q^p'^-^ t)^L {u, q, r) {t - rf-' + . 
The fractional stationary condition ([5]) reduces to 

BT 

P'^'^ (t) = - ^ {u (r) , g (r) , r) {t - rr~' ; (9) 
the first equation in ([7]) to 

u(r)=i?"-''q(T) ; (10) 
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while the second equation in ([7]) takes the form 

uq 

Using (0) and (HH]) in Ull) we arrive to: 
'dL 



— {u{r),q{T),r) (t - rf-^ 



(11) 



D 



/3,Q 



9m 



(12) 



Simple calculations show that is equivalent to ([3]), that is, Theorem 13. II is 
a generalization of Theorem 12.21 to the fractional optimal control problem ^ . 

Remark 3.2. Let us define the Poisson bracket of two dynamical quantities / 
and g, with respect to coordinates q and fractional-momenta p"''^, by 



{/,5} = 



df dg df dg 



dp"'f^ dq dq dp°'''^ ' 
The fractional Hamiltonian system ([7]) can be written in the following form: 

D^^^qir) = {H'^^^q} , 
D^J^.^.p'-'^ (t) = {H'^'^p^^P} . 

Proof. Theorem 13.11 is proved applying the (a, /3) fractional Euler-Lagrange 
equations to the augmented action integral ((S]), i.e. applying ^ to 



Via) 



where H"''^ — H"'^ (u (t) , q {t) ,p"''^ (t) The Euler-Lagrange equation 

with respect to q gives the second equation of the fractional Hamiltonian system 
([7]); the Euler-Lagrange equation with respect to u gives the fractional station- 
ary condition ([5]); and, finally, the Euler-Lagrange equation with respect to p"'^ 
gives the first equation of ([7]). □ 

In classical mechanics, constants of motion are derived from the first integrals 
of the Euler-Lagrange equations. In the fractional case it is necessary to change 
the definition of constant of motion in a proper way ^1 1101 [TT] . Here, in order 
to account the presence of the fractional decaying friction force F"':^ ([4]), we 
propose the following definition of fractional constant of motion. 

Definition 3.1. We say that a function C of r is a fractional constant of motion 
if and only if D^i'^.^C = 0. 
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Corollary 3.2. If L and do not depend on q, then it follows from the frac- 
tional Hamiltonian system (cf. Theorem \3 . 1\) that D^'^.^p"'^ (r) — 0, i.e., p"^'^ 
is a fractional constant of motion. 

Definition 3.2 (fractional momentum of order (a,/?)). Associated witti an 
{a, (3) fractional action-like variational problem ^ we define the fractional mo- 
mentum by ©-(Uni): 



Corollary 3.3. For the fractional problem of the calculus of variations ^ 
with §1=0, the fractional momentum of order is a fractional constant 

of motion. 

4 Conclusions 

We generalize previous results of the fractional variational calculus by using 
the action-like variational approach together with Ricmann-Liouville fractional 
derivatives of order (a, /?). We claim that the (a, (3) fractional action-like varia- 
tional problem here introduced offers a better mathematical model for weak dis- 
sipative and nonconservative dynamical systems. Both Lagrangian and Hamilto- 
nian approaches are considered. We derive the (a, (3) fractional Euler-Lagrange, 
Hamilton and Poisson equations. Standard Noetherian constants of motion are 
violated due to the presence of a fractional decaying friction term. To solve 
the problem we introduce a new notion of fractional constant of motion. The 
class of fractional Hamiltonian systems thus obtained has two parameters and 
is wider than the standard class of Hamiltonian dynamical systems. Our ap- 
proach is different from the ones found in the literature: different from S. Muslih 
and D. Baleanu approach [14 , where no fractional integral action is used and 
no decaying friction term considered; different from the recent approach [20] of 
A. Stanislavsky, where the generalization of the classical mechanics with frac- 
tional derivatives is based on the time-clock randomization of momenta and the 
coordinates are taken from the conventional phase space. 
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